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Well-posedness and exact controllability of the mass balance 
equations for an extrusion process 

Mamadou DIAGNE* Peipei SHANG 1 ^ and Zhiqiang WANG* 


Abstract 

In this paper, we study the well-posedness and exact controllability of a physical model 
for a food extrusion process in the isothermal case. The model expresses the mass bal¬ 
ance in the extruder chamber and consists of a hyperbolic Partial Differential Equation 
(PDE) and a nonlinear Ordinary Differential Equation (ODE) whose dynamics describes 
the evolution of a moving interface. By suitable change of coordinates and fixed point 
arguments, we prove the existence, uniqueness and regularity of the solution, and finally 
the exact controllability of the coupled system. 
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1 Introduction 

The analysis of moving boundary problems has been an active subject in the last decades 
and their mathematical understanding continues to be an important interdisciplinary topic 
for various engineering applications. Representive complex physical systems describing bio¬ 
logical phenomena and reaction diffusion processes such as Stefan problems in cristal growth 
processes are still calling for many open questions related to their exact controllability and 
the design of highly efficient output feedback control laws for stabilization purposes. Among 
these challenging problems, one can mention the swelling nanocapsules studied in [2j, the 
lyophilization process applied to pharmaceutical industry mm, the cooking processes de¬ 
scribing the volume change in food material IMl, the mixing systems (model of torus reactor 
including a well-mixed zone and a transport zone) and the Diesel Oxidation Catalyst (DOC) 
presented in ['29J. 
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In general, the key point in such problems is to find a suitable change of coordinates which 
transforms the system with moving interface into a system defined on a fixed domain. While 
leading to more complicated infinite dimensional operator, such approaches allow to utilize the 
existing results that establish the existence, uniqueness, regularity and continuous dependence 
of solutions for many systems of conservation laws. For the well-posedness problems, we refer 
to |T1 l3l fl5l fl6l 21, 28j (and the references therein) in the content of weak solutions to systems 
(including scalar case) of conservation laws, and [24112T ] in the content of classical solutions 
to general quasi-linear hyperbolic systems. The control problems for hyperbolic conservation 
laws have been widely studied for a long time. For controllability of linear hyperbolic systems, 
one can see the important survey m- The controllability of nonlinear hyperbolic equations 
(or systems) are studied in (5j EJ EJ [23)[26j]- Moreover, provides a comprehensive survey 
of controllability and stabilization in PDEs that also includes nonlinear conservation laws. 

In this paper we consider the well-posedness and exact controllability of the Cauchy 
problem for a physical model of the extrusion process which describes the mass transport 
phenomena in an isothermal extruder chamber. Mathematically, the process is described by 
a hyperbolic PDE defined on a time-varying domain. The dynamics of the spatial domain 
is governed by an ODE expressing the conservation of mass in the extruder and the PDE 
expresses the convection phenomenon due to the rotating screw. More detailed description 
of the model is given in Section [2j We mention that the first result concerning the mathe¬ 
matical analysis of the extrusion model as transport equations coupled via complementary 
time varying domains is proposed in m, where the well-posedness for the linearized model 
of the extruder is obtained by using perturbation theory on the linear operator. 


Our proof of the well-posedness (see Theorem 3.1) of the Cauchy problem for the ex¬ 
trusion model relies on the characteristic method and the fixed point arguments. The ir¬ 


regularity of the solution is proved as well (see Theorem 3.2), which is useful when one 


considers the asymptotic stabilization of the corresponding closed-loop system with feedback 
controls m- In this context the stabilization of a non-isothermal food extrusion process 
including the temperature and the moisture content dynamics is investigated in [14] . The 


idea to prove Theorem 3.3 is to construct a solution to (2.2)-(2.4) which also satisfies the 


final conditions. The way of such construction is based on the controllability result of the 
linearized system together with fixed point arguments (see for example 0) 

The organization of this paper is as follows: First in Section [2j we give a description 
of the extrusion process model which is derived from conservation laws. The main results 
(Theorems [371} [372} [373} concerning the well-posedness, regularity and the controllability 
of the normalized system are presented in Section [3} while their proofs are given in Sections 
[4} [5} [6] respectively. 
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2 Description of the extrusion process model 


Extruders are designed to process highly viscous materials. They are mainly used in the 
chemical industries for polymer processing as well as in the food industries. An extruder is 
made of a barrel, the temperature of which is regulated. One or two Archimedean screws are 
rotating inside the barrel. The extruder is equipped with a die where the material comes out 
of the process (see Fig. [I]). 

feed moving heater 



0 PFZ l(t) FFZ L 


Fig. 1: Schematic description of an extruder 


In an extruder, the net flow at the die exit is mainly due to the flow of the material in the 
screw axis direction. The die resistance influences highly the transport along the extruder and 
induces an accumulation phenomenon towards the down barrel direction. The accumulation 
of the material permits us to represent the mass balance in the extruder using the filling 
ratio along its screw channel. More precisely, the spatial domain of the extruder might be 
partitioned into a Fully Filled Zone (FFZ) and a Partially Filled Zone (PFZ). The flow 
in the FFZ depends on the pressure gradient which appears in this region due to the die 
restriction. The PFZ corresponds to a conveying region that is submitted to the constant 
atmospheric pressure and the transport velocity of the material depends only on the screw 
speed and pitch. These two zones are coupled by an interface which is located at the spatial 
coordinate where the pressure gradient changes from zero to a nonzero value. Basically, the 
moving interface evolves as a function of the difference between the feed and die rates. In 
the sequel, the spatial domain of the extruder will be taken as the real interval [0, L\, where 
L > 0 is the length of the extruder. Let us denote by fit) G [0, L\ the position of the thin 
interface, the domain of the PFZ is then [0, fit)] and the FFZ is defined on [fit), L] (see Fig. 

0 - 

Considering the following change of variables m 


r4i/ = —— in PFZ and x e->- y 
y fit) y 


x — fit) 
L - fit) 


in FFZ 


( 2 . 1 ) 


respectively, the time varying domains [0 ,l(t)] and [l(t),L\ can be transformed to the fixed 
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domain [0,1] in space. For the sake of simplicity, we still denote by x the space variable 
instead of y. More precisely, we consider the problem for the corresponding normalized 
system defined on (0, T) x (0,1). 

Defining the filling ratio along the PFZ spatial domain, namely, f p (t,x) as a dynamical 
variable [20|[22j, then the mass balance in this region is written as follows: 


dtf P (t , x) + a p d x f p (t, x ) = 0, in M+ x (0,1), 

< f P (0,x) = fp(x), in (0,1), (2.2) 

Jp{t, 0) = p a VeffN(t) ’ in M + , 


where 


a p (x,N(t),l(t),f p (t,l)) = 


CN(t)-xF(l(t),N(t),f p (t,l)) 

l(t) 


(2.3) 


and Op is the transport velocity of the material, ( is the screw pitch, F tn denotes the feed rate, 
po is the melt density, V e jf is the effective volume and N(t) is the rotation speed of the screw. 
F is the dynamics of the moving interface described by equations (2.4)-(2.5). The interface 
motion is generated by the gradient of pressure which appears in the FFZ [18, T3J [20] , 22j 
and under the assumption of constant viscosity along the extruder (the isothermal case), its 
evolution is given by the following 


where 


with 


i(t) = F(l(t),N(t),f p (t, 1)), inK+, 

m = i° , 

F(l{t), N (f), f p (t, 1)) = N(t)g(l(t), f p (t , 1)), 


9W)Jp{^ !)) = 


_ C f P (t, i) 

[B Po + K d (L - l(t))} (1 - f p (t, 1)) 1 - f p (t, 1)' 


c K d (L-l(t)) 


(2.4) 


(2.5) 


( 2 . 6 ) 


In (2.6), K d denotes the die conductance and B is the geometric parameter. 

In the whole paper, unless otherwise specified, we always assume that 1° £ (0, L), f p £ 
VF 1,oo (0,1), Fi n ,N £ L°°(0, T), Fi n /N £ hF 1 ’ oo (0,T). For the sake of simplicity, we denote 
from now on ||/||l<x= (\\f\\ w i,oo, ||/|| L 2, resp.) as the L°° (IF 1 ’ 00 , L 2 , resp.) norm of the 
function / with respect to its variables. 


3 Main Results 


In this section, we present the main results on the well-posedness and exact controllability 


of the coupled system (2.2)-(2.4), we have the following two theorems. 


Theorem 3.1. Let T > 0 and (l e , N e , f pe ) be a constant equilibrium, i.e., 

F(l e ,N e ,fpe) = 0 


(3.1) 
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with 0 < f pe < 1, 0 <l e < L. Assume that the compatibility condition at (0,0) holds 

Fin ( 0 ) _ f0/n\ 

p o V eff N(0) 

Then, there exists eo > 0 (depending on T) such that for any e E (0, £o], if 


\l°—le\ + ||/p(0 — fpeWw 1 ’ 00 + 


Fin(-) 


p 0 V eff N(- 


~ fpeWw 1 ’ 00 + ||-W(') — Ne\\L°° 


(3.2) 


(3.3) 


Cauchy problem (2.2)-(2.4) admits a unique solution (l, f p ) E VF 1,oo (0, T)xW 1,oo ((0, T)x(0,1)), 
and the following estimates hold 


IK(') — leWw 1 ’™ < C £ o • £, 
ll/p('j ') — fpeWw 1 -™ — Fe 0 ■ £, 

where C £0 is a constant depending on £q, but independent of e. 


(3.4) 

(3.5) 


Theorem 3.2. Under the assumptions of Theorem 3.1, we assume furthermore that /))(•) E 

F ( 

H 2 { 0,1), — —— E H 2 (0,T), and the compatibility condition at (0,0) holds 


PoV e ffN(- 


(fO) (o) + ■ — ( _ Fin ^ _^ 

Up)xU <ZN(0) dt \poV eff N(t)) 


t =o 


= 0. 


(3.6) 


(3.7) 


Then, there exists eo > 0 (depending on T) such that for any e E (0, eo], if 

\l°~le\ + \\fp{-)~ fpeWmpdp) + II p Q V™fN(-) ~ fp e W H2 (.°’ T ) + I! FI (') -^e 11 L°° 

Cauchy problem (2.2)-(2.4) has a unique solution ( l,f p ) E W 1,oo (0,T) x C°([0, T]; H 2 {f), 1)) 
with the additional estimate 


II/p(v) ~ fpe\\c°([o,T]-,H 2 (o,i)) < C £Q • e, 
where C £0 is a constant depending on eo, but independent of e. 


(3.8) 


Remark 3.1. The solution in Theorem \3.1\ or in Theorem 3.2 is often called semi-global 
solution since it exists on any preassigned time interval [0,T] provided that ( l,f p ) has some 
kind of smallness (depending on T), see 


Remark 3.2. We have the hidden regularity that f p E C°([0,1]; H 2 (0, T)) in Theorem 
(see QEf for the idea of proof). 


3.2 


The problem of exact controllability for Cauchy problem (2.2)-(2.4) can be described as 
follows: For any given initial data (Z°, fp(x)), any final data (l 1 , f 2 (x)), to find a time T and 
controls Fi n (t ) and N[t) such that the solution to Cauchy problem (2.2)-(2.4) satisfies 

l{T) = l\ (3.9) 

f P (T,x) = /». (3.10) 

Our result is the following theorem on local controllability in the sense that the initial and 


final data are both close to the given equilibrium determined by (3.1). 
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Theorem 3.3. Let 


T P : = 


l. 


(N e 


(3.11) 


be the critical control time. Then, for any T > T e , there exists v\ > 0 suitably small such 
that, for any v G (0, v\\, l°,l 1 G (0, L) and fp,fp € VF 1,oo (0,1) with 

\l° — le\ + 1 1 1 — le I + II fp(') ~ fpeWw 1 ’ 00 + II fp(') ~ fpeWw 1 ■°° < (3.12) 

there exist N G L°°(0,T) and Fi n G L°°(0,T) satisfying 

Fin(-) 


p 0 V e ffN(- 


~ /pellvK 1 ' 00 + 11^(0 — ||Z/°° < 0,1 • I/, 


(3.13) 


such that the weak solution (l(t), f p (t,x)) to Cauchy problem (2.2)-(2.4) satisfies the final 
condition (3.9)-(3.10). Here, C Ul is a constant depending on v\, but independent of v. 


4 Proof of Theorem 3.1 


In order to conclude Theorem 3.1, it suffices to prove the following lemma on local well- 


posedness of Cauchy problem (2.2)-(2.4). 


Lemma 4.1. There exist £\ > 0 and 6 > 0 suitably small, such that for any e G (0,£i], 
1° G (0, L), /° G W 1,oo (0, 1), F in /N G W 1 - oo (0,T) with 


\l°~le\ + \\fp(-) — fpe\\w 1 ’°°F 


Fin(-) 


PoVeffN(- 


~ fpeWw 1 ’ 00 F\\N(-) — N e \\L°o<£, 


(4.1) 


Cauchy problem (2.2)-(2.4) admits a unique local solution on [0, <5], which satisfies the follow¬ 
ing estimates 


\l{t ) - l e \ < C £ i ■£, Vt G [0, <5], 

II fp(t,-) - fpe IIvr 1 ’ 00 < C e 1 ■£, Vt G [0, J], 


(4.2) 

(4.3) 


where C £1 is a constant depending on £i, but independent of £. 


Let us first show how to conclude Theorem 3.1 from Lemma |4.1[ By Lemma |4.1| we 
take £2 e (0, £i] such that C £l ■ £2 < £1 • Then for any £ G (0, £ 2 ] and any initial-boundary 


data such that (4.1) holds, Cauchy problem (2.2)-(2.4) admits a unique solution on [0,<5]. 
Furthermore, one has 


I l(S) - l e | < C £1 ■ £ < £ 1 , 
||/p(^> ') /pellw 1 ’ 00 — Cfii ‘ ^ — ^1- 

By taking (1(5), f p (6, •)) as new initial data and applying Lemma 
of Cauchy problem 


4.1 


(4.4) 

(4.5) 

on [5, 25} , the solution 


-(2.4) is extended to [0,25]. For fixed T > 0, we can extend the 


local solution to Cauchy problem (2.2)-(2.4) to [0, T] eventually by reducing the value of £q 
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and applying Lemma 4.1 in finite times (at most [T/5\ + 1 times). Therefore, to conclude 


Theorem 3.1, it remains to prove Lemma |4.1| 

Proof of Lemma |4.1f The proof is divided into 4 steps. 

Step 1. Existence and uniqueness of (/(•), f p (-, 1)) by fixed point argument. 

Let Si > 0 be such that 


□ 


Denote 


0 < Si < min{Z e , L l e , _/p e , 1 fpe}- 


|E||vyi,oo : = ^ esssup \D a F(xi, X 2 , 0 : 3 )|, 

I I 0<X1<L 

N e -ei<X2<N e +ei 

0<a:3<l 

= 1)), t G [0,T], 


(4.6) 


(4.7) 


(4.8) 


For any given 5 > 0 small enough (to be chosen later), we define a domain candidate as 
a closed subset of C°([0, <5]) with respect to C° norm: 

0 5 , ei := {d/ e C°([0,5])|4r(0) = (Z°,/°( 1)), ||T(-) - (l e , f pe )\\ c o {m ) < d}- (4.9) 

We denote by £(s;t,x), with (s,£(s;t,x)) G [0, t\ x [0,1], the characteristic curve passing 
through the point (t, x) G [0,(5] X [0,1] (see Fig. [ 2 ]), i.e., 

x ^ ^ 10 ^ 

f(t: t, x) = x. 

Let us define a map 5 := (SibS^), where £ : fZ, 5 )£l —> C°([0, (5]), *F 1 —> 5 r (4') as 

Um) --=1°+ f F(l(s),N(s)J p (s,l))ds, ( 4 . 11 ) 

Jo 


Mm) := /p(£(0; t, 1 )). 


(4.12) 


Solving the linear ODE (4.10) with a p given by (2.3), one easily gets for all <5 small and all 
0 < s < t < 5 that 

ft f(i(g),Af(g-),/p(g-,i)) d<j ft (N(a) r HWM fAsA)) ds 

£(s;t, 1) = e Js l( - a ) — / • e Js da. (4.13) 

Js H 0 ") 

It is obvious that 5 maps into £l$ !ei itself if 

le S 1 le £l L l e £1 


0 < 6 < min 


in {T, 


}• 


/-(AT I \ ’ II p|| 1 || p|| [ j (4-14) 

Q(i\l e -rS 1) \\r llw 1 -’ 00 ||L Ww 1 ^ ' 

Now we prove that, if 5 is small enough, 5 is a contraction mapping on f 2 < 5 £1 with respect 

to the C° norm. Let = (/, / p ), T = ( l , f p ) G fi,s j£l . We denote by £(s; f, x) the corresponding 

characteristic curve passing through (t, x): 

= ap(£(s-,t,x),N(s)J(s),f p (s, 1)), ^ 1 


£(t; t, x) = x. 
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(4.16) 


Similarly as (4.13), one has for all 5 small and all 0 < s < t < 5 that 


ft F(l( a ),N(a)Jp(<T,l)) ■ (N(a) r F d<.s),N( S )J p (s,l)) d 

£(s;i, 1) = e Js 1 . • e Js da. 


1(a) 


Therefore it holds for all t E [0,5] that 


|3i(*)(*) - = f F{I(s),N(s), f p (s, 1)) ds - f F(l(s), N (s), f p (s, 1)) ds 

Jo Jo 


^^IlFlI^i.ooUT - ^||cO([0,5])- 


(4.17) 


On the other hand, it follows from (4.12), (4.13) and (4.16), that for all t E [0,5], 

\um)-um)\ 


<»/; 


px 


I L°° 


ft F(l(<r),N(<T)J p (tr,l)) ^ , t F(l(a),N(<r),f p (a,l)) , 

*10 l(c r) _ gdO Z(<t) 


+ 


(N(a ) C^(cr) f ; ds 

/» ' C * *(<r) ' 6 


dcr). 


By (4.7) and the fact that 4 / ,T E 0<5. £] and l(t) > l e — ei > 0, |lV(f)| < iV e + £i, it follows 
that for all t E [0,5], 


i um) - u*m\ < cs \\$(•) - fpeiiw^w* - ^Hcom, 


(4.18) 


where C > 0 is a constant independent of (\h, x I'). Finally, combining (4.17) and (4.18), we 
can choose 5 small enough such that 


1„ = 


II^W - t?(^)||c 0 ([0,<5]) < 9 11^ - ^lie 0 ([0,5])■ 


(4.19) 


Banach fixed point theorem implies the existence of the unique fixed point (Z(-), f p (-, 1)) of 
the mapping T = 5(\h) hr ^<5,ei- 

Step 2. Construction of a solution by characteristic method. 

With the existence of (/(•), / p (-, 1)) and 5 > 0 by Step 1, we can construct a solution 


to Cauchy problem (|2.2|)-( 2.4). For every (t,x) in [0,5] X [0,1], we still denote by £(s;i, x), 


with (s, £(s; t,x)) E [0, t] x [0,1], the characteristic curve passing through the point (t, x), see 


(4.10). Since the velocity function a p is positive, the characteristic £(s;t,x) intersects the 
x-axis at point (0 ,/3(t,x)) with f3(t,x ) = £(0;i, x) if 0 < £(f; 0,0) < x < 1; the characteristic 
£(s;t,x) intersects the t- axis at point (r(t, x),0) with g(r(t, x);t, x) = 0 if 0 < x < £(t; 0,0) 
(see Fig. [2]). Moreover, we have (see m Lemma 3.2 and its proof, Page 90-91] for a more 
general situation) 


dr(t, x) 
dx 

dP(t,x) 

dx 


= -Kr(t,x)) 
CN(r(t,x)) 

ft F(l(a-),N(a-) : f p (tT,l)) , 
— gdQ Z(cr) 


(4.20) 

(4.21) 
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t 



o p(u) 


Fig. 2: The characteristics £(s;t,x) and f3(t,x), r(t,x). 


We define f p by 

f P (t,x) = < 


W if 0 < a: < 0,0) < l,0<t<6, 

poV e ffN(T{t,x)) 

fp(f3(t, x)), if 0 < £(t; 0,0) < x < 1, 0 < t < 5. 


(4.22) 


Then it is easy to check that (l,f p ) £ IT 1,oo (0, 5) xPT 1,oo ((0, 5) x (0,1)) under the compatibility 


condition (3.2) and (l, f p ) is indeed a solution to Cauchy problem ( |2.2| )-(2.4). 

Step 3. Uniqueness of the solution. 


Assume that Cauchy problem (2.2)-(2.4) has two solutions (l, f p ), (l, f p ) on [0, <5] x [0,1]. it 
follows that (Z(-), f p (-, 1)) = (Z(-), f p (-, 1)) since they are both the fixed point of the mapping 
5: T = 3 r ('k) in f2«5 i£l . This fact implies that the characteristics x) and x) coincide 
with each other and therefore so do the solutions f p and f p by characteristic method. 

Step 4. A priori estimate on the local solution. 


By definition of f p and assumption (4.1), it is clear that for all t £ [0,<5], 


II ~ fpe\\L°° < £■ 

Thanks to 

|/(f)| =\F (l(t),N(t ), f p (t , 1)) - F (l e , N e , fpe) I 

<\\F\\ w l,oo(\l(t) - l e | + \N(t) - N e | + | f p {t, 1) - f pe |). 
<\\F\\\yl,oo\l(t) — Z e | + 2e||F||^yl.oo, 


(2.4), (O), (4.23) and assumption (4.1), we get for all t £ [0, <5] that 


(4.23) 
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which yields (4.2) from (4.1) and Gronwall’s inequality. On the other hand, from (4.22) 
dfp 


dx 


< 


< 


d / F in (T(t,x)) 
dx \p 0 V e ffN{T(t,x)) 
F in (-) 


+ 




PoV e ffN(- 




pe 


w 1 


dr 


dx 


+ 


/£(■)-/, 


pe 


w 1 


d/3 


dx 


(4.24) 


Combining ( 4.20[ ), (4.21), (4.24) and assumption (4.1), we obtain (4.3) which concludes the 
proof of Lemma 4.1 □ 


5 Proof of Theorem 3.2 


Before proving Theorem 3.2 let us recall a classical result on Cauchy problem of the 
following general linear transport equation 


ut + aft, x)u x = b(t, x)u + c(t, x), (t, x) 6 (0, T) x (0,1), 

u(0,x) = uo(x), re G (0,1), 
u(t,0) = h(t), t€(0,T), 

where aft, x) > 0, a, a x , b £ L°°((0, T ) x (0,1)) and c € L 2 ((0, T ) x (0,1)). 


(5.1) 


We recall from [T, Section 2.1], the definition of a weak solution to Cauchy problem (5.1). 


Definition 5.1. Let T > 0, uq € L 2 ( 0, 1), h G L 2 (0,T) be given. A weak solution of Cauchy 


problem (5.1) is a function u £ C°([0, T]\ L 2 (0,1)) such that for every r G [0,T], every test 
function (p G C' 1 ([0, T] x [0,1]) such that ip(t, 1) = 0, \/t G [0, T\, one has 

rr rl , . rl 


u[d t p + ad x p + (a x + b)p\ + cp^J dx dt + / u(t, •) dx 


Jo Jo 

— f uo<p(0,-)dx— f ha(-, 0)<p(-, 0) dt = 0. 

Jo Jo 


(5.2) 


We have the following lemma 

Lemma 5.1. Let T > 0, uq G L 2 (0,1) and h G L 2 (0,T) be given. Then, Cauchy problem 


(5.1) has a unique weak solution u G C°([0, T]; L 2 (0,1)) and the following estimate holds: 

ll' u llc°([0,T];L 2 (0,l)) < C'(II' u o||l2 (0i i) + ||^||l 2 (0,T) + II c IIl 2 ((0,T)x(0,1)))> ( 5 - 3 ) 

where C = C(T, ||a||i,°o((o,T)x(o,i)); ll a a:||L°°((o,T)x(o,i))) PIU oo ((o,r)x(o,i))) a constant inde¬ 
pendent of uq, h, c. 


For the proof of Lemma 5.1, one can refer to m for classical solution or m Theorem 
23.1.2, Page 387] for Cauchy problem on M without boundary. 

Proof of Theorem |3.2[ By Theorem |3.1| and Lemma |5.1| it suffices to prove that the 
systems of f Pxx satisfies all the assumptions of Lemma 5.1 
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Differentiating (2.2) with respect to x once and twice give us successively that 


dtf Px (t,x) + a p {t,x)d x f Px (t,x) = ~a Px {t,x)f Px (t,x), in (0,T) x (0,1), 

— in (0,1), 

-l(t) d / F in (t) 


fp*(0,x) = f° x (x) 


fpM, °) = 


(N(t) dt \ poV ef fN(t)) ’ 


and 


in (0, T), 


(5.4) 


•W 0,x) = f$ xx (x), 
fpxAt, 0 ) = 


C N(t) L 


x )dxf Pxx {t,x) = -2 a Px 

(t, 

v)f Pxa (t,x), 

in (0, T) x (0,1), 

), 



in (0,1), 

rF(l(t),N(t),f p (t,l)) 

d 

( F in (t) \ 

(5.5) 

I CN(t ) 

dt 

\p 0 V eff N(t)J 



d ( l(t) d ( F in (t ) 


dt\(N(t) dt\p 0 V eff N(t) 


in (0, T), 


with 


^ (N(t)-xF(l(t),N(t),f p (t,l)) u , —F(l(t), N(t), f p (t, 1 )) 

a p {t,x) = -—-, a Px {t,x) = - 


m 

From the assumptions that f p E FI 2 (0,1), 


m 


F- 

J- 1.7 


(5.6) 


PoVeffN(-) 


E H 2 (0,T) and the compatibility 


conditions (3.2) and (3.6), one easily concludes Theorem 3.2 by applying Lemma 5.1 


Cauchy problem (5.5) 


to 

□ 


6 Proof of Theorem 3.3 


The idea to prove Theorem 3.3 is to construct a solution to Cauchy problem (2.2)-(2.4) 


which also satisfies the final conditions. The way of such construction is based on the con¬ 
trollability result of the linearized system together with fixed point arguments (see 0) 

For any fixed initial data (l°,f p (x)) and final data (l 1 , f p (x)) close to the equilibrium 
(lafpe), we define a domain as a closed subset of C 0 ([0,T]) with respect to C° norm: 

:= {$ e C°([0, T])|<f>(0) = (1°, fp(l)), <F(T) = (l 1 ,/}(!)), ||<&(-) - (l e , f P e)\\c°([o,T]) < £i}, 


where the constant £\ > 0 is determined by (4.6). We study the exact controllability for the 


linearized system deduced from equations (2.2)-(2.4), replacing (l(t),f p (t, 1)) by (a(t),b(t)) 


in functions F and a p : 


i(t) = F a ’\t ), 

dtf P (t , x) + a. p b {t, x)d x f p (t, x) = 0, 
= f P (0,x) = f p (x), 


f P (t, 0) = 


*Z b 0t) 


Po V eff N^(t)' 


in (0, T), 

in (0, T) x (0,1), 

in (0,1), 

in (0, T), 


( 6 . 1 ) 
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where 


F a ’\t) = N a ’ b {t)g a ' b (t), 
a,b, t s = C R d (L - a(t)) 


C*p b (t,x ) = 


_ (b(t) 

[Bp 0 + K d (L - a(Z))] (1 - b(t)) 1 - b(t) ’’ 

(N a ’ b (t ) - xF a ’ b (t ) 
a(Z) 


( 6 . 2 ) 

(6.3) 

(6.4) 


By assumption T > T e , it is possible to find controls N a,b (t ) and F°^ b (t) such that the solution 
of (6.1) satisfies (3.9)-( phT0 ). 

For any = (a, b) £ f l Sl,T , we choose the control function N a,b (t ) as 


/i _/0 1 

N a,b (t) := ——- Z£[0,T], 

W T 5 a ’ 6 (t) J 


such that 




thus Z(Z) is a linear function of Z: 


Z 1 -/ 0 


Z(t) = Z u H-^Z £ [0, T]. 


By assumption (3.12), it follows then 


\l(t) - Z e | = 1(1 - |)(Z° - l e ) + f -{l l - Z e )| < 2^ < 2^. 


(6.5) 


( 6 . 6 ) 


(6.7) 


( 6 . 8 ) 


Next let us construct the desired control F ? “' b (t). For this purpose, we shall give the 
expression of f p (t,x) by characteristic method. We denote by ^ a,b (s;t,x) the characteristic 
passing through (t,x) £ [0,T] X [0,1]: 


U^ b (s-t,x) _ b b ^ _ (N a ’ b (s)-^ b (s-t,x)F“’ b (s) 


ds 


= a“’ (s,£ a ’ (s\t,x)) = 


a(s) 


£ a,b (Z; t, x) = x. 


(6.9) 


Suppose that the characteristic £ a,b (s ; 0, 0) intersects the line x = 1 at (Zq , 1). Then for every 
Z £ [0,Zg’ b ], the characteristic passing through (t, 1) intersects with x-axis at (0,/3 a,fc (Z)); for 
every t £ (Zq’ 6 , T], the characteristic passing through (Z, 1) intersects with Z-axis at (r a,6 (Z), 0). 
The characteristic passing through (T, 1) intersects with Z-axis at (Z“’ b , 0) (see Fig. [ 3 ]). Hence, 
we have 


e a ’ 6 (Zg’ 6 ;0,0) = 1, £ a,6 (0; Z, 1) = /3 a,b (t), C’ b {r a ' b (t)-t, 1) = 0, £ a ’ b (Z“’ b ; T, 1) = 0. 
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► X 


Fig. 3: The characteristics £ a,6 (s;i, x) and Zq’ & , (3 a,b (t), T a,b (t), f“’ 6 . 


Solving (6.9) and using (6.4) and (6.6), we obtain 


J 0 Ta(a) aa _ 


<’ b (N a ’ b (a) r^ds r 

- - — ■ e Jo Ta (°) dcr = 1, 


Jo a(a) 

ft Z 1 —Z® ft C ((j} fCT Z 1 —Z® j n , 

fo Ta(rr) d(7 _ 1 -. Jo Ta(s) dS faj = fi a ’ b (t) , 

Jo a{a) 

(N a ' b (cr) r, l ^ds 
/ -——• e J - a ’ b M T ”W da = 1, 

J T a ’ b (t) a \ a ) 


t-t P-i 0 A _ rt 

e J-r a ’ b {t) Ta(a) a(J _ 


f t °-, b Ta.(v) da f 1 (N a ' ( a ) fa,b Ta(s) , 

eh — / - , ; -eh da = 1. 


ds 


t a a b a ( fJ ) 


We define the control function F°^(t) as the following 

F“’ & ( t ) = { /l(t)p ° K//iVQ ’ 6(0, tE[0X ^ 

\f^ a ’ b (T;t,0))poV eff N a ’ b (t), t G [tf,T], 

where h(t) is any artificial W 1,a ° function satisfying 

I !HJ — fpeWw 1 ’™ < v l 
and the following compatibility conditions 


/i (0) = /p(0), h(t°f b ) = fl( 1) 


( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 


(6.14) 


(6.15) 

(6.16) 
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Inspired by (4.22), we define a map T := (.Fi,Fh) : fi £l,T —> C°([0,T]), 4> (->• F(<f>) as 


T 1 ^)(t):=l(t) = l 0 + l ^J-t 1 t E [0, T], 
l/i(T a > 6 (t)),tE(# 6 ,T], 


(6.17) 

(6.18) 


where fg’ b , /3 a,b (t) and T a,b (t ) are defined by (6.10), (6.11) and (6.12) respectively. 

Now we prove that T is a contraction mapping on f2 £l,T provided that v\ > 0 is small. 
Obviously, the existence of the fixed point implies the existence of the desired control to the 
original nonlinear controllability problem. 

Let 1/1 < ei/2. By ([6781), ([6171) and (Ogb, F($) E fi ei > T for all 4> E n ei - T . 


For any 4> = (a, b) E Q ei,T , we denote by ^ a,b (s',t, x) the characteristic passing through 


(t,x) E [0, T] X [0,1] as in (|6.9[) upon replacing a, b by a, b. Correspondingly, we define t\ 


a,b 
0 > 


j3 a,b (t), T a,b (t ), t“’° as in (|6.10[), (|6.11[), (|6.12[) and (|6.13[) upon replacing a, b by a, b. 


,a,b 


By definition of T = [T\, F 2 ) in (6.17) and (6.18), we get J-\ ( ( I>) = Fi($) and thus 


||^(*F) -^(4>)|| c o ([0 ,t ]) = ||F 2 ($) -F 2 (<b)|| c o ([ o, T]) = sup |(F 2 ($) - F 2 ($))(t)|. 

te[o,T] 

Without loss of generality, we may assume that t^ b > t^ b . Hence, we need to estimate the 
point-wisely |F 2 ($)(t) — F 2 ($)(i)| on the time interval [0, t^}, [tQ ,b ,t^’ b ], [tQ ,b ,T] respectively. 
For any given t G [0,tg’ b ], by (3.12) and (6.18), we have 


1(^(4) - Jb(*))(f)l = - fp{P a,b (t)) 

<II&IM/ 3 a ' S (t) -/ 3 “'*(t) 
<n\ /r f (t) - 


By (6.11), it is easy to get 

I P a,b (t) - f3 a,b (t )| < C{\\a — a||c°([o,x]) + ||£> - &||c°([o,T])) < C*||4> - $||co([ 0 ,t])> 

Here and hereafter, we denote by C various constants that are independent of ($, <F). There¬ 
fore, 


sup |(F 2 (T) - F 2 (<f>))(f)| < Cv i||4> - < I , || C '0([o,T])- 


(6.19) 


tG[0,*S’ 

,a,b 


For any given t E [iq , T], by (6.15) and (6.18), we obtain 

mm - Kmm = 1 h(r~ a ’\t)) - h(r a ’ b m < ^(t) -r^wi. 


( 6 . 20 ) 


By flop , for every t E [t^ ,T], 

£ a,6 (s; t, 1) = ^ d<T ~ J s 


d , l CN“’ b (a) jyl=£ ds 


a{a) 


. gJ s Ta(s) s dcr, s E [r a,b (t),T]. 
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Therefore 


sup |£ a,b (s; t, 1) — £ a,b (s; t, 1)| < C(\\a - a|| c o [0iT] + ||&-&||c°[o,T])- 

s£[T a ’ b (t),T] 


Then, for any t E [fo’ b ,T], we have 


T a - 6 (t) - T°- e (t) | < : 


a.b ( 


inf oi 


a,b 


rr a,b (t) 

/ a a /(s;C' b (s;t,l))ds 


inf ct 


a,b 


r a ’ b (t) 


- a a /(s-C’ b (s-t,l))) ds 


a,b( pa,b( 


where 


<< 7 ( 1(0 - a||c , 0([0,T]) + \\b - ft||c°([o,T]) 

<(711$ - $||co([0,T])> 


inf a“’ b := inf a“’ 6 (f, x) > 0 

p (t,x)e[o,T]x[o,i] p 


(a,b) 6f! e l’ 


is independent of $. Therefore, 


sup |(J r 2 ($) - J r 2 ($))(t)| < Cv i||$ - $||c°([0,T])- 


te[4’ b ,T] 


For any given t E [tg’ fe ,£o’ 6 ], by (6.16) and (6.18), we have 


1 j- 2 ($) - J 2 ($)i =1 m^ b {t)) - Hr a ’ b m 


ip' 
rO ( oa,b / 


<1 W“'“(0) - J?(0)| + \h(T‘-‘(t)) - A(0)| 
<ll/?Jli-K S {<)l + ll'i<lk»|T“’ 1 (t)l 
< I - 1 (|/3 a ’ S (t)| + |r“-‘(i)|). 


Note that for any t E [^q’ & , ^q’ 6 ] , 


_ t- r* d-i° r 1 dN a ’ b (rr) r° l l -l° j„ 

C’ (s; t, 1) =e fs j • e f ‘ T da, s E [0, t]. 

Js a(a) 

tra,b/ n . + 5 ,,b s _ 


Then, using £ a,6 (0; £q , 1) = 0, it follows that 


r >6 coi = ir ,6 (o; t, 1) - T’ 6 (o ; c°, i)i < c\t - c°i < c\%° - 


ca,b( n . .a,b 


a.b 1 


,a,b ,a,6i 


( 6 . 21 ) 


( 6 . 22 ) 


(6.23) 


(6.24) 


(6.25) 


On the other hand, by (6.12) and note that T a ’ b (t^' b ) = 0, we can prove that for any t E 

r,a ,6 ,a, 6 -i 

1*0 ’ r o J 5 

(6.26) 


|r a,6 (t)| = | T a ' b {t) - T a ’ b (t a 0 ' b ) | < (7|f - to’ b | < C|£“’ b - i?’ 6 ' 


By definition of £q and £q’ , similarly to the derivation of ( |6.21| ) and (6.22), we get 


\C - C°l < (7(||a - a||c°([o,T]) + \\b - i>||c°([o,T])) < C'll^ _ ^||c°([o,t])' (6.27) 
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Then it follows from (6.24), (6.25), (6.26) and (6.27), that 


sup |(^ r 2 ( < h) - J r 2 (^’))(t)| < Cis i||<h - 4 > ||co([o,t])- 
te[t a 0 ’ b ,tf] 


Combining (6.19), (6.23) and (6.28), we can choose v\ small enough such that 


ll-^W - •^ 7 ( < ^ > ) llc°([0,T]) = ll-^^) - ^2(^)|bo([0,T]) < 2 11^ “ < ^ > llc°([0,T])- 


(6.28) 


(6.29) 


Then the contraction mapping theorem implies that T has a unique fixed point (Z(-), f p (-, 1)) 
in n £l,T . Therefore we find a solution (Z, f p ) to Cauchy problem (2.2)-(2.4) which also satisfies 
the final conditions ( |3.9[ )-( 3.10). Moreover, because of the uniqueness of solution, the desired 
control function N(t) and Fi n (t ) can be chosen by substituting the fixed point (Z(-),/ p (-, 1)) 
into (6.5) and (6.14) respectively. This concludes the proof of Theorem 3.3 □ 
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